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Abstract 


The simulation of radio frequency interference (RFI) cancellation by applying a spatial filtering technique for 
phased array feed (PAF) is presented. In order to better reflect the characteristics of PAF, a new signal model is to 
add the coupling coefficient among elements of PAF to the conventional array signal model. Then the subspace 
projection (SP) algorithm is used to cancel RFI from the correlation matrix of the signal, and finally, the 2D power 
image is drawn. The power variation of signal-of-interest direction and RFI direction before and after using the SP 
algorithm is analyzed. The new signal model and simulation strategy can be used to test and verify the beamformer. 
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1. Introduction 


A phased array feed (PAF) receiver, as an example of the 
next generation microwave receiving system, has been widely 
studied all over the world. However, the applications of spatial 
interference cancellation using PAF are rare in radio astronomy 
(Wang et al. 2021). At present, most radio frequency 
interference (RFI) mitigation methods are for narrow-band 
and pulsed or intermittent interfering signals, these are based on 
a time-frequency character of the signal and the RFI are 
normally blanked. In this paper, the simulation of RFI 
mitigation based on spatial information using the PAF receiver 
is presented. A signal model with planar coupling effect is 
established, and the effect on the main beam is analyzed. The 
subspace projection (SP) algorithm is used to simulate PAF 
receiver cancellation RFI and the performance of output power 
is studied in signal-of-interest (SOD and RFI directions. 

The SP algorithm is used to simulate PAF mitigation RFI 
and performance of output power in different directions is 
studied. Different from conventional array, PAF is a tightly 
distributed array designed for reflector antennas. The conven- 
tional array signal model cannot include a coupling effect 
among elements of PAF. This paper presents the performance 
of the PAF coupling effect in data correlation operation from 
baseband data. Then, radio telescope pattern and RFI 
cancellation are simulated and analyzed by combining antenna 
simulation software, the MaxSNR and SP algorithms. 

The Mk II PAF receiver developed for ASKA by the 
Commonwealth Science and Industrial Research Organization 
in Australia has been compressed from 0.8—1.74GHz to 
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1.2-1.75 GHz due to interference from mobile phone and 
digital television signals after Parkes and Effelsberg tests 
(Chippendale et al. 2016). The American FLAG project 
considered adaptive filtering technology to mitigate RFI in 
beamforming and verified the performance of Linear Constraint 
Minimal Variance (LCMV) and SP algorithms in PAF 
(Warnick et al. 2007; Landon et al. 2012). Now, there is no 
PAF receiver fully put into use in China, but several research 
projects and groups, including the Five-hundred-meter Aper- 
ture Spherical radio Telescope (FAST) team of the National 
Astronomical Observatories (Chai et al. 2020) and Qi Tai 
Telescope (QTT) project (Wang 2014), are rapidly carrying out 
research. 

QTT, the largest fully steerable single dish radio telescope 
currently under construction in the world, has been selected by 
Nature as the most noteworthy science project in 2023.5 Its 
prime focus will be equipped with a 96-element dual- 
polarization PAF array. (Ma et al. 2019), where beamforming 
will be implemented in GPU of digital terminal (Pei et al. 
2022). There is no doubt that this PAF receiver is the most 
complex system in QTT, and its development is full of 
challenges. However, a suitable numerical model can be used 
for theoretical analysis, validation, or improvement of 
beamforming algorithms, etc. Therefore, a baseband numerical 
model and algorithm simulation that can reflect the character- 
istics of the PAF receiver become meaningful. Ivashina et al. 
(2011) optimized the joint field of view and beam rotational 
symmetry of PAF receiver through numerical simulations in 
the Westerbork Synthesis Radio Telescope (WSRT) project. 


> https: //www.nature.com /articles/d41586-022-04444-3 


Research in Astronomy and Astrophysics, 24:015019 (8pp), 2024 January 


Deep-space 


Interfering [ - 
N object 


satellite [] 


Spatial 
null 


Multiple 
beams 


PAF signal 


High gain primary antenna 


Figure 1. The main advantages of PAF are that multiple beams can be formed 
simultaneously to provide joint field of view, actively control the number and 
shape of the beams, and spatial cancellation of interference. 


Landon et al. (2012) studied the improvement of the SP 
algorithm for moving RFI in numerical simulation. The 
coupling effect of PAF was not considered in these signal 
models. Nybo Jeffrey (2019) established a signal model 
containing coupling characteristic in the FLAG project for 
the test and verification of GPU beamforming, but the signal 
model is a simple linear model, which is not suitable for planar 
array. 

With the use and development of radio in people’s 
production and life, there is an urgent need to deal with the 
ubiquitous interference, such as television signals, FM radio 
transmission, air navigation equipment. Even radio-quiet zones 
are not immune to man-made radio interference, such as 
satellite downlinks. As shown in Figures 1 and 2, PAF in radio 
astronomy provides a promising method. The spatial features 
of RFI can be extracted by sampling the focal plane field to 
track and remove interference without discarding data. Even 
interference sources that completely cover the astronomical 
signal spectrum can be canceled. This technique is expected to 
improve the capability of radio telescopes in the presence of 
interference. 

The remainder of this paper is organized as follows. 
Section 2 provides the new model of PAF, detailed definition, 
and mathematical preparation of the SP beamforming. In 
Section 3 we customize a simulation strategy and verify the 
performance of new signal model in the main beam on the 
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QTT-PAF model. The SP algorithm is used to cancel the data 
containing RFI in the beamformer, and power variation before 
and after using the SP algorithm is analyzed. Section 4 is a 
summary. 


2. Preliminaries 
2.1. Signal Model of PAF in Time 


As depicted in Figure 2, assuming narrowband operation of 
an M element PAF, the M x 1 complex baseband data vector at 
time n sample is given as 


x[n] = c(as[n] + vd[n] + n[n]) (1) 


where s[n] is the SOI with array response vector a (normalized 
spatial characteristics of celestial bodies in PAF, sometimes 
called array manifolds, including effects of the reflector dish 
and the geometry of PAF), d[n] is the single interfering source 
with corresponding array response v (usually not in the 
direction of main focus). In general, a is constant because 
telescopes are mechanically controlled at the same line of sight 
as SOI. As the interferer is random or the dish tracks SOI, v is 
changing with time. c is the M x M matrix, which describes the 
coupling characteristic of PAF, and the following section will 
introduce in detail. Noise vector m[n] is assumed to be spatially 
non-white, Gaussian, and the correlation arises mainly from the 
PAF coupling effect. Signals s[n], d[n], and n[n] are assumed 
to be temporally wide sense stationary. Beamformer output y is 
the weight w = [w], W2,....Wyl” sum of signals received by 
each array element and is computed as 


y= w?x. (2) 


2.2. Signal Correlation Analysis 


Usually, the signal sampling rate is much faster than the 
movement of interference, so v is approximated to a constant in 
a short-time integration (STI). The SOI and RFI are supposed 
to be uncorrelated from each other. These should also be 
independent of the noise of each receiving chain. The 
covariance matrix of the signal in the jth STI window of 
length T samples is estimated to be 

: 1 DT-1 1 P 

R, =— x[n]x" [n] = —X;X; 

T > [n]x" [n] A 
= œ (aRa” + vR” + N) Ea (3) 
Couple 


Rsor Rrry Rwoise Couple 


where: 


1. X = [xT], xT + 1),--xl(i DT 1]] is the 
M x T matrix, the voltage output of PAF through the T 
sample points observation. 

2. N is the M x M noise covariance matrix, supposed to be 
diagonal. We assume N = XI, where X% is noise power 
and IJ is the identity matrix. 


Research in Astronomy and Astrophysics, 24:015019 (8pp), 2024 January 


Fronted 


Wu, Li, & Chen 


Backed 
beamforming 


X2 


reflector 
couple 


7 RI 


noise 


Da XM AE 


Figure 2. Generalized representation of PAF receiver. 


3. R, is the 1 x 1 SOI covariance matrix, which is also the 
power of SOI in an STI window. The generated SOI has 
Gaussian, centered, white and stationary characteristics. 

4. R, is the Q x Q RFI covariance matrix. Q is the number 
of interferer. 

5. c is the M x M matrix, which is the coupling coefficient 
among the elements. The tight arrangement among 
elements of PAF increases the radiation bandwidth and 
destroys the beam pattern. The coupling phenomenon 
means that when a wave arrives at PAF, the induced 
current in one element will radiate some energy back to 
other elements in PAF. There is more induction among 
the near elements and less induction among the far 
elements. In the model in Figure 3, the elements are 
placed in the rectangular coordinate system, and the 
distance between two elements can be expressed by a 
tuple (Ax, Ay). At this time, when the signal of PAF is 
correlated, the coupling effect caused by the distribution 
of the elements can be assumed to be a covariance matrix 
€, which is expressed as follows: 


in 2 2 2 
Oo {^ 0,0 T AVio gai Ký AY 


o Attot Ary gaita 


gyArin+t Aym 
TE 2 
Oo Ariat AY (4) 


o | Axi m+ Aym 


where ø is the coefficient of the coupling effect (generally 
determined by the inherent characteristics of the elements 
material). So c is the Cholesky decomposition matrix 
of €. 


€= 


o Aot Ayo o Ahat Ay. 
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Figure 3. The layout of PAF elements in the plane. 
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2.3. The SP Algorithm 


The mathematical principle of SP projection in the PAF 
receiver can be summarized as a linear transformation from one 
vector space to another. 

Compared with the waveguide feed, the PAF receiver can 
extract the spatial characteristic (the position information) of a 
target signal, and the signal can be selectively mitigated or 
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Figure 4. PAF model equipped with the QTT project consists of 96 x 2 
elements. The selected elements in the figure are used for one beam. 


accepted, which provides a new idea to cancel RFI for the 
single-aperture radio telescope. In this section, we will explain 
the working principle of the SP algorithm from the perspective 
of vector space of data. In an STI window, the expected value 
of the signal is estimated by the covariance matrix, and the 
spatial feature vector of RFI is extracted by the array 
information processing method. § is used to represent the 
output signal space of PAF. According to different radio 
astronomical signal processing targets, S can be divided into 
three subspaces, namely SOI subspace, RFI subspace, and 
noise subspace. The three subspaces are expressed as follows: 


Ra) = O, 
R([vi,..., vol) = V, 
R(n) = N = R(x[n]) = S. (5) 


where R (-) represents the column space. We want to project the 
output signal of PAF into a special space W, which satisfies that 
W is perpendicular to V and W belongs to S. The matrix 
orthogonal projection method can be used, the expression is as 
follows: 


P- = I — vw(v8y)|y4 (6) 
which is such that P-v = 0. Therefore, when the spatial filter is 
applied to the data covariance matrix, all the energy due to 


interference will be nulled. The column space of P* can be 
described as 


R(P+) =W=02-V (7) 


where the column space of J is the whole space Q = R°. It can 
be understood as finding a set of vector basis that does not 
contain RFI spatial feature, so that the projection value of some 


Wu, Li, & Chen 


GRASP-E 


Jb *grd file 


Python 


= o 


Figure 5. Flow chart of simulation strategy. GRASP-E simulation results are 
read by Python in the data format of *.grd. 


data containing RFI spatial feature on this set of vector basis is 
Zeros. 

Then, the P+ transform is performed on PAF weight w to 
obtain the weight vector w; with RFI mitigation capability. 
The expression is as follows: 


Plw = w. (8) 


In the linear transformation of Equation (8), P+ vertical has an 
important feature: idempotency, which determines that the 
projection only plays a role in the first mapping and does not 
superimpose multiple mapping effects. 


3. Simulation 


In this section, the simulation design is based on the PAF 
model in QTT. As shown in Figure 4, a detailed numerical 
model is used for 32 of the 96 elements. This scheme refers to 
the design of Ma et al. (2021). In this numerical model, we 
generate the data model described in Equation (1) and use that 
to realize the beamforming simulation and performance 
analysis of the SP algorithm canceling RFI. 


3.1. Simulation Conditions 


The tools we use in this simulation are the GRASP-E and 
Python. The former is a free version of a general purpose 
reflector antenna designed for students by TICRA company in 
Denmark,° and the latter is an open source data processing and 
drawing script language. The whole simulation strategy can be 
divided into four parts, and the flow chart of the simulation 
strategy is shown in Figure 5. 

Main Reflector. GRASP-E is used to simulate the pattern of 
110 m reflector. As the specific design details are not disclosed, 
the influence of the support structure is not considered in this 
simulation. Figures 6 and 7 list the reflector model and pattern. 


© https: //www.ticra.com /software / grasp-student-edition / 
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Figure 6. QTT 110 m paraboloid model, focal diameter ratio f/d = 0.33. The 
feed is conical and taper angle is 74°29. 


Signal Generator: The complex exponential signal of PAF 
including SOI, RFI, and noise is generated (Li et al. 2023), 
which satisfies the description in Equation (3) statistically. The 
main parameters of the Signal Generator in this simulation 
process are (a) the element’s position parameter table (which 
largely determines the shape of the beam); (b) the proportional 
coefficient between the element’s spacing and operating 
frequency; (c) the ratio of sampling rate to signal frequency; 
(d) the location information of SOI and RFI in space; (e) the 
SNR of SOI and RFI, where SNR = 10log,, >, P, is signal 


n 
power and P,, is noise power. 


SP-Beamformer: The main purpose of this part is to calculate 
the vector P% and to realize beamforming. P% is the weight of 
MaxSNR Wymaxsnr reweighted, the expression is w; = 
Pwmaxsner. It mainly includes signal covariance calculation, 
matrix decomposition, and other operations. The important 
parameters considered in this simulation process are (a) the size 
of STI window; (b) the normalized spatial characteristic of RFI. 
The error of measured estimated ¥ and v plays a key role in the 
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Figure 7. Normalized far-field pattern, resolution 55 x 55. The half-power 
beamwidth (HPBW) of QTT is 0°025 at 1.5 GHz. 


cancellation of RFI. (v is given directly for convenience, and 
this paper does not involve the measurement and estimation 
of v). 

Visualization of Power: Radio astronomy often focuses on 
the measurement of celestial radiation power. In an STI of 
PAF, the average power can be measured by digital beam 
control. The expression is as follows: 


P = wRw (9) 


where w is the normalized spatial direction vector. The power 
measurement in an STI can be realized by using this formula to 
traverse the imaging area. 


3.2. The Influence of Coupling Matrix in PAF 
Beamforming 


In order to analyze the influence of coupling coefficient in 
PAF beamforming, 32 of 96 QTT-PAF elements are selected to 
beamforming one beam. The distribution of these elements is 
shown in the selected ones in Figure 4, which are located in the 
center of the PAF and have good directional characteristic. The 
Signal Generator is used to generate an observation signal with 
a power far lower than that of noise and without RFI, which is 
to represent the scene when the telescope deep space 
exploration in a non-interference environment. The detailed 
parameters of the Signal Generator are shown in Table 1. 

The Signal Generator generates two kinds of data, one with 
a coupling coefficient and the other without. As shown in 
Figure 8, the correlation matrices for these data are plotted 
separately, along with the correlation matrix for the coupling 
coefficient. 
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Figure 8. The correlation matrices of the two kinds of data generated by the Signal Generator, where STI = 4000. (a) Correlation matrix without a coupling 
coefficient in the data; (b) Correlation matrix when ø = 0.4; (c) Correlation matrix when the coupling coefficient is included in the data. 


MaxSNR ki couple 


-0.45 -0.22 


$/(0.01*) 
¢/(0.01*) 


0.0 


1.0 1.0 1.0 
ae So 
0.45 -0.22 0.22 0.45 -0.22 (0) 
0.45 0.45 0.45 
0.8 0.8 0.8 
0.22 0.22 
0.6 - 0.6 - 0.6 
(0) d 0 d 
0.4 0.4 0.4 
-0.22 -0.22 
0.2 0.2 0.2 
-0.45 0.45 


0.01") 


0.0 0.0 


Figure 9. There are the 2D power images of the telescope in the data condition of Table 1. (a) shows the power image without the coupling matrix in the data. (b) The 
power image with coupling proof in the data set is displayed. (c) shows the power difference between the first two images. 


Table 1 
SOI Simulation Parameters 


Parameters SNR(in)/dB SNR(out)/dB 6/° o/° f/GHz 
SOI —30 —30.99 0 0 1.5 
RFI null null null null null 
noise white noise 


From the correlation results in Figure 8, it can be seen that 
the coupling coefficient enhances correlation among elements 
while retaining the autocorrelation of conventional array. We 
use MaxSNR beamformer to complete beamforming of the 
results shown in Figure 8. In the ideal matching condition, the 
telescope pattern is obtained by weighting the Main Reflector 
and PAF pattern. We selected the area around the main lobe to 
draw 2D images using Visualization of Power, as shown in 
Figure 9. 


It can be seen from Figure 9, when compared with data 
images without a coupling coefficient, the main lobe shape of 
the data images with a coupling coefficient is deformed, 
rotation symmetry is lost, and the first “null” depth is 
shallower. This is consistent with our understanding of the 
mutual coupling effects of PAF. 


3.3. Subspace Projection Cancellation RFI 


Subspace projection is an algorithm to control the “null” 
position of a telescope pattern. In the presence of interference, 
the SP algorithm uses projection operator P to cancel out RFI in 
the correlation matrix R, estimated by M x M dimension in 
an STI. 

Normally in radio astronomy observation, the main lobe of a 
telescope is aligned with the target source, and the interference 
source enters through the sidelobe. The parameters in Table 2 
are used to simulate the PAF array signal with coupling effect 
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Figure 10. The SNR of the beamforming signals in the SOI and RFI directions. Here, the purple lines are the result of the MaxSNR beamformer; the red lines are the 
result of the SP-Beamformer; the cyan lines are the output of the MaxSNR beamformer with pure SOI (SOI and noise). 


Table 2 
SOI and RFI Simulation Parameters 
Parameters SNR(in)/dB = SNR(out)/dB 0/° o/°  f/GHz 
SOI —30 —29.86 0 0 1.5 
RFI 20 20.67 —0.6 0.6 1:5 
noise white noise 


by the Signal Generator when interference is present. It is 
assumed that there is only one interference source and its 
direction is random (excluding SOI direction). 

We use the SP-Beamformer and MaxSNR beamformer to 
achieve beamforming against the elements signal, respectively. 
Then, the relationship between SNR and STI is calculated 
respectively in the direction of SOI and RFI, and the results are 
shown in Figure 10. In the results, we found that the presence 
of RFI in the beamforming has an effect on the power 
performance in the direction of SOI, and this effect is related to 
the size of STI. This is different from what we understood at 
the beginning. 

In Figure 10(a), compared with the pure SOI results, the 
presence of RFI causes the energy to climb in the SOI 
direction. The main reason for this is the leakage of energy 
from the RFI direction to the SOI direction. After the 
beamforming of SP-Beamformer, the energy climb is 
improved, and this ability tends to a fixed value with the 
increase of STI size. In Figure 10(b), RFI is canceled by “null” 
after SP-Beamformer processing. The result of this deletion has 
a large dynamic range and will not be improved with the 
increase of STI size, but it can meet the requirement of RFI 


cancellation. In general, the SP algorithm can well meet the 
cancellation requirements of RFI, but it cannot fully restore the 
power of SOI, which is caused by va = 0. 

In conclusion, an appropriate increase of STI size can make 
the performance of the SP algorithm to restore SOI more stable, 
although the algorithm cannot be completely restored. For 
environments with strong RFI, the size of STI has less effect on 
cancellation because RFI dominates the signal power. There- 
fore, in the design of beamformer, the STI size should consider 
the stability of power restore in the direction of SOI while 
taking into account the performance of GPU. 

Figure 11 is the sky power image of QTT, implemented by 
Visualization of Power. In Figure 11(a), the image after RFI 
failure can be shown that the RFI from the sidelobe is the main 
response area, while SOI response at the central position is 
covered up. Figure 11(b) is the image recovered by SP 
technology. It can be seen that RFI is completely eliminated 
and the response of SOI is highlighted. Figure 11(c) is the 
difference between the two figures. 


4. Conclusion 


The main content of this paper are as follows: (1) From the 
perspective of data correlation, the array data model which can 
reflect the coupling characteristic of planar PAF is established, 
and the coupling effect can lead to a loss of rotational 
symmetry and a shallower first “null”. (2) The simulation of the 
combination of PAF and reflector antenna pattern is completed. 
The SP algorithm is used to cancel the RFI in PAF receiver and 
draw 2D power response images. Calculate the change of SOI 
power in RFI direction before and after using the SP algorithm. 
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Figure 11. There is the 2D power image of QTT in the data condition of Table 2, STI = 2048. 
Although the simulation of the model mentioned in the paper is References 


mainly focused on the digital back end, without considering 
LNA, antenna efficiency, and other factors, the data model and 
simulation strategy can be used in the test to verify the 
beamformer and shorten the development time of the PAF 
receiver. 
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